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> ■ ABSTRACT 

I We present new evidence for mature stellar populations with ages > 100 Myr in mas- 

OO . sive galaxies (M s t c iiar > 10 10 M Q ) seen at a time when the Universe was less than 

CO 1 1 Gyr old. We analyse the prominent detections of two z w 6 star-forming galaxies 

(SBM03#1&#3) made at wavelengths corresponding to the rest-frame optical using 
the IRAC camera onboard the Spitzer Space Telescope. We had previously identified 
these galaxies in HST/ ACS GOODS images of Chandra Deep Field South through the 
"z'-drop" Lyman break technique, and subsequently confirmed spectroscopically with 
^ ■ the Keck telescope. The new Spitzer photometry reveals significant Balmcr/4000A 

Oh" discontinuities, indicative of dominant stellar populations with ages > 100 Myr. Fit- 

q . ting a range of population synthesis models (for normal initial mass functions) to the 

' HST /Spitzer photometry yields ages of 250 — 650 Myr and implied formation redshifts 
£0 j Zf ~ 7.5 — 13.5 in presently-accepted world models. Remarkably, our sources have 

d • best-fit stellar masses of 1.3 — 3.8 x 10 10 M Q (95% confidence) assuming a Salpeter 

IMF. This indicates that at least some galaxies with stellar masses > 20% of those 
of a present-day L* galaxy had already assembled within the first Gyr after the Big 
1 Bang. We also deduce that the past average star formation rate must be comparable 

5-j \ to the current observed rate (SFRtjv ~ 5 — 30Af Q yr _1 ), suggesting that there may 

have been more vigorous episodes of star formation in such systems at higher redshifts. 
Although a small sample, limited primarily by Spitzer 's detection efficiency, our result 
lends support to the hypothesis advocated in our earlier analyses of the Ultra Deep 
Field and GOODS HST/ ACS data. The presence of established systems at z as 6 
suggests long-lived sources at earlier epochs (z > 7) played a key role in reionizing the 
Universe. 
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1 INTRODUCTION 

In recent years, advances in detector efficiency, large ground- 
based telescopes and space-based observatories such as the 
Hubble Space Telescope (HST) and Spitzer Space Telescope, 
have revolutionized studies of the high redshift Universe. 
Searches based on Lyman-a emission are now at last uncov- 
ering many galaxies at 2 < z < 7 (e.g., Kodaira et al. 2003; 
Rhoads et al. 2004). The Lyman-break technique (Steidel, 
Pettini & Hamilton 1995; Steidel et al. 1996; 1999) has like- 



wise proved successful in selecting high-redshift star-forming 
galaxies. This utilises the rest-UV continuum break seen 
shortwards of Lyman-Q and caused by H 1 absorption in the 
intergalactic medium. The redshift range z m 6 is of great 
importance, as it heralds the end of the reionization of the 
Universe (Becker et al. 2001; Kogut et al. 2003) which might 
be achieved through star formation. Using HST and the new 
Advanced Camera for Surveys (ACS; Ford et al. 2003), the 
Lyman-break technique has been pushed to this early epoch 
(Stanway, Bunker & McMahon 2003; Bouwens et al. 2004a; 
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Yan & Windhorst 2004; Giavalisco et al. 2004) by using the 
i' and z' filters to isolate i'-drop galaxies. Ground-based 
follow-up spectroscopy (Bunker et al. 2003; Stanway et al. 
2004a,b; Dickinson et al. 2004) has shown that a colour cut 
of (i' — z')ab > 1.5 mag reliably finds star-forming galaxies 
at z ~ 6 with modest foreground contamination (primar- 
ily from low-mass stars, and passively evolving galaxies at 
z ~ 1 — 2). Subsequent low-dispersion slitless spectroscopy 
with HST/ACS (Malhotra et al. 2005) confirms the nature 
of i'-drops as z ~ 6 galaxies with a small contaminant frac- 
tion of stars and low-redshift sources. 

Yet an important part of the puzzle is missing: the 
i'-band drops are selected in the rest-frame UV, and are 
therefore known to be actively forming stars. However, it 
is unclear whether these objects suffer from significant red- 
dening due to dust (in which case the star formation rates 
will have been underestimated), or if there is an underly- 
ing older stellar population which has been recently rejuve- 
nated. Publicly-available Spitzer imaging with the Infrared 
Array Camera (IRAC; Fazio et al. 2004) as part of the Great 
Observatories Origins Deep Survey 1 (GOODS; Dickinson & 
Giavalisco 2003; Dickinson et al. m prep) allows us to ad- 
dress both questions. 

A key benefit of Spitzer photometry arises because the 
IRAC camera samples wavelengths longwards of the age- 
sensitive Balmer & 4000A breaks at z ~ 6. Accordingly, 
we have analysed the IRAC images of GOODS-South field 
(centred on the Chandra Deep Field South, Giacconi et al. 
2002), where we have already selected i'-drop galaxies from 
the GOODS-South ACS images (Stanway, Bunker & McMa- 
hon 2003) and the Ultra-Deep Field (UDF; Bunker et al. 
2004). 

Our goal in this paper is thus to focus on the in- 
frared properties of four z ~ 6 galaxies for which there 
are robust spectroscopic redshifts, based on Lyman-a emis- 
sion. The properties of the entire i'-drop population in 
the GOODS fields is considered in a forthcoming paper 
(Eyles et al. in prep). Our four confirmed sources include 
two sources from Stanway, Bunker & McMahon (2003): the 
brightest confirmed i'-drop in the GOODS-South field (with 
z' — 24.7 mag) SBM03#3 (with a spectroscopic redshift of 
z = 5.78 from Keck/DEIMOS, Bunker et al. 2003); and the 
brightest i'-drop in the UDF, SBM03#1 with z = 25.3 mag 
(spectroscopically confirmed with Keck/DEIMOS by Stan- 
way et al. 2004a as having a redshift of z = 5.83). Both of 
these spectroscopic redshifts were independently confirmed 
by Dickinson et al. (2004). The other two sources come 
from the Gemini Lyman-a at Reionization Era (GLARE, 
Stanway et al. 2004b) spectroscopy with Gemini/GMOS: 
GLARE#3001 (z = 5.79, z' = 26.4 mag) and GLARE#3011 
(z = 5.94, z = 27.2 mag). 

A plan of the paper follows. In Section|2]we describe the 
Spitzer imaging data and the methods used to fit stellar pop- 
ulations to the broad band photometry derived collectively 
from Spitzer, HST ACS & NICMOS images, and ground- 
based near-infrared data. We discuss the implications of the 
age and stellar mass estimates in Section [2] Our conclusions 
are presented in Section |3] Throughout we adopt the stan- 



see http : //www. stsci . edu/f tp/science/goods/ 



dard "concordance" cosmology of Am = 0.3, J1a = 0.7, 
and use Ho = 70 kms -1 Mpc -1 - in this cosmology, the 
Universe today is 13.67 Gyr old, and at z = 5.8 the age 
was 992 Myr. All magnitudes are on the AB system (Oke & 
Gunn 1983). 



2 OBSERVATIONS, DATA REDUCTION AND 
STELLAR POPULATION FITTING 

2.1 Spitzer Observations 

We concentrate here on the shortest-wavelength Spitzer im- 
ages of GOODS-South taken with IRAC as part of the "Su- 
per Deep" Legacy programme (PID 194, Dickinson et al. in 
prep). We have analysed the first half of the GOODS-South 
Spitzer data, taken with 39 Astronomical Observation Re- 
quests (AORs) observed between 8-16 February 2004. 

The IRAC camera comprises four channels, each with 
a 256 2 InSb array imaging a 5.2' x 5.2' field with a pixel 
size of ~ l'.'22. Images were taken through four broad-band 
infrared filters, with central wavelengths at approximately 
Acent = 3.6 /im, 4.5 /im, 5.6/^m and 8.0 /im (channels 1-4), 
and widths of AAfwhm = 0.68, 0.87, 1.25, 2.53 /im respec- 
tively. The individual frame times for each channel were 
200 s (except those taken at 8.0 fim, which comprise 4 inte- 
grations of 50s at each position). Over the course of the 39 
AORs, a large 2x2 mosaic of pointings was executed, with 
smaller random sub-dithers when the pattern was repeated, 
giving a 10' x 10' coverage for each filter. Each AOR com- 
prised 10-11 pointings. Channels 1 & 3 (3.6 (im & 5.6 fim) 
have the same pointing, offset by 6.7' from the common 
pointing of channels 2 & 4 (4.5 jum & 8.0 pm). As a result, 
only a portion of the field is observed in all four wavebands. 
Once data from the second epoch (when the telescope has 
rotated by 180 degrees) has been gathered, images of the 
entire 10' x 16.5' field matching the GOODS-South ACS 
survey (Giavalisco 2003) in all four filters will be available. 
Meanwhile, by design, the overlap region common to all fil- 
ters covers the Ultra Deep Field (Beckwith et al. 2003). The 
total exposure time in each channel is « 86 ksec, depending 
on location. 



2.2 Data Reduction 

For our investigation, we used the pipeline-processed IRAC 
images at the 'Post-Basic Calibrated Data' (PBCD) stage 3 , 
details of which can be found in the Infrared Array Camera 
Data Handbook, Version 1.0 (Reach et al. 2004). The main 
steps in the pipeline include dark current subtraction, ap- 
plication of flatfields, flux calibration (in units of MJy/sr) 
and mosaicing of the individual frames of each AOR after 
application of a distortion correction. PBCD sets have a re- 
fined pointing solution to an accuracy of 0'.'2 derived from 

2 Adopting the first-year WMAP best-fit cosmology from Spergcl 
et al. (2003) with Q M = 0.27, Q A = 0.73, and H = 
71 kms -1 Mpc -1 makes negligible difference (the masses and lu- 
minosities would be 4 per cent greater, and the Universe would 
be 4 per cent younger at z = 5.8). 

3 The PBCD images are available from 

http : //ssc . spitzer . caltech . edu/popular /goods/ irac_ 04/ 
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2MASS point source catalogue objects in the field of view. 
Outlier rejection is performed during the mosaicing process. 

We used IRAF to combine the 39 PBCD mosaics (one 
for each AOR), which were taken at different telescope roll 
angles. The IRAF task wregister was used to rotate the 
frames to a common roll angle. The individual registered 
frames were scaled to the same integration time and all the 
frames were offset by the median counts in the centre of 
the area (which was common to all frames) to account for 
floating bias levels, and the 39 frames were combined using 
imcombine (weighting by the exposure times). During this, 
residual unrejected cosmic rays, bad pixels, and the column 
pull-down and muxbleed detector effects were removed us- 
ing a percentile clipping method "pclip", rejecting at the 
3cr level. The PBCD images were converted from units of 
surface brightness (MJy/sr) into flux units of fiJy per pixel 
by multiplying the data units by 10 12 x (pixel solid angle) : 
a numerical factor of 34.98. The images were matched to 
the vl.O reduced z'-band tiles of the GOODS-South field 
released by the GOODS team 4 . The astrometry was found 
to be consistent to within ~ 0'.'2. While our reduction was 
underway, the GOODS team released an enhanced dataset 
(DR1) 5 , employing a 'multidrizzle' technique similar to that 
used successfully on HST/ ACS GOODS data. This provides 
combined images with a pixel scale of « 0'.'6. The magni- 
tudes listed in this paper are determined from this "driz- 
zled" data, and we have used our independent reduction 
as a consistency check. Upon cross-checking the astrome- 
try, photometry and noise properties, they were found to be 
consistent to 0.1 mag (for bright sources) and 0'.'2. 

In the final GOODS-South co-added DR1 'drizzled' im- 
ages, we measured the FWHM of the PSF to be w l'/5 
in channels 1 & 2 (3.6 & 4.5 /jm), and « l'/8 in channel 3 
(5.6 /im) & ~ 2'.'l in channel 4 (8/xm). Even for the short- 
wavelength IRAC images, the Spitzer PSF is much larger 
than the typical size of the z ~ 6 i'-band drop galaxies 
(whose half-light radii r M < 0"2, Bunker et al. 2004). At 
Spitzer resolution, these galaxies are clearly unresolved and 
so we treat them as point sources. The galaxy images in the 
various wavebands are shown in Figures [5] El & ID 

To construct spectral energy distributions of our 4 
spectroscopically-confirmed i'-band drop galaxies in the 
GOODS-South field (detailed in Section and Table [TJ we 
undertook aperture photometry in the various images. In 
order to maximize the signal-to-noise ratio (S/N) and min- 
imize possible confusion with other foreground objects, we 
used a diameter « 1.5 x FWHM for the IRAC images, appro- 
priate for unresolved objects. The aperture diameters were 
4, 4, 5 & 6 'drizzled' pixels for the 4 channels (3.6, 4.5, 5.6 & 
8.0 /im), corresponding to 2'.'4, 2'.'4, 3'.'0, & 3'.'7 respectively. 
We used the IRAF digiphot . phot package to measure the 
enclosed flux at the exact coordinates determined from the 
GOODSvl.O and UDF i'-band images, taking the residual 
background from an annulus between 12" and 24" radius. 
We applied aperture corrections to compensate for the flux 



4 Available from 

ftp: //archive . stsci . edu/pub/hlsp/goods/vl 
a For the enhanced dataset GOODS DR1 see 
http : //data. spitzer . caltech. edu/popular/goods 
20050209_enhanced_vl/Documents/goods_drl .html. 
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Figure 1. The fractional flux enclosed within a circular aperture 
as a function of aperture diameter, determined from 4 bright but 
unsaturated isolated point sources in the IRAC channel 1 (3.6 £tm) 
GOODS-South image. The small aperture used for faint sources 
(2'.'4-diameter) is marked: this encircles 52 per cent of the energy, 
implying an aperture correction of 0.7 mags. Note that the official 
Spitzer calibration for IRAC uses a 24"-diamctcr aperture, so the 
fraction of the flux there is set to 1.0 by definition. 



falling outside the aperture: these were w 0.7 mag, as deter- 
mined from bright but unsaturated point sources in the im- 
ages, with apertures of diameter 18" in channels 1 & 2, 24" 
in channel 3, and 30" in channel 4. The curve of growth for 
four stars is shown in Figure Q (note that the official Spitzer 
calibration for IRAC also uses a similar 24"-diameter aper- 
ture, so the aperture correction there is 1.0 by definition). 
These aperture corrections are consistent with those derived 
for the First Look Survey (Lacy et al. 2005). 

Checks were made to ensure the objects were not con- 
taminated by the presence of any neighbouring bright fore- 
ground sources. An object was deemed to be uncontami- 
nated if no bright source lay within a 4" radius from it. All 
four objects satisfied this criteria; the bright source close to 
SBM03#3 lies on the 4" limit (Figure^. The Spitzer magni- 
tudes in Table|5|have the neighbour subtracted out (i.e., are 
deblended) but this has only a ~ 5% effect on the measured 
flux in our 2'.'4-diameter aperture; less than the error bar. 
The subtraction of the neighbour was done using the "GAL- 
FIT" software (Peng et al. 2002) . This takes into account the 
PSF (we give it a PSF for the drizzled image derived from 
stacking 4 bright but isolated stars) and a best-fit galaxy 
surface brightness profile. An exponential disk worked best 
in the case of this neighbour, but it is not very sensitive to 
the exact functional form with the poor resolution of Spitzer. 
Figure 0] shows the postage-stamp image of SBM03#3 with 
the model fit of the neighbour subtracted. 

The noise for each of the four channels was checked in 
two different ways. First, we derived an estimate based on 
a Poisson model using the detector gain, number of frames 
combined, and the background counts (adding back the zo- 
diacal background estimate subtracted by the pipeline but 
recorded in the header) . Secondly, we measured the standard 
deviation in background counts of the images. As the mo- 
saicing process introduces correlations between pixels, we 
also made noise estimates using the individual BCDs and 
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ID RA (J2000) Dec (J2000) Redshift Age of Universe / Gyr 



SBM03#1+ 03:32:40.01 -27:48:15.0 5.83 0.95 

SBM03#3 03:32:25.61 -27:55:48.7 5.78 0.96 

GLARE#3001 03:32:46.04 -27:49:29.7 5.79 0.96 

GLARE#3011 03:32:43.10 -27:45:17.6 5.94 0.93 

t this is also known as GLARE#1042 (Stanway et al. 2004b), #20104 (Bunker et al. 2004), SiD002 (Dickinson et al. 2004) and YWOla 

(Yan & Windhorst 2004). 

Table 1. The four spectroscopically-confirmed i'-drop galaxies in the GOODS-South Field. 



assuming it decreased as the square root of the number of 
frames. These estimates lead to 3<r limiting AB magnitudes 
of 26.5 and 26.1 in 2'.'4 apertures in channels 1 and 2, re- 
spectively, and 23.8 and 23.5 in 3'.'0 apertures in channels 
3 and 4, respectively. There will be additional background 
fluctuations caused by faint galaxies (i.e. confusion noise), 
which will increase the noise. Both methods produce consis- 
tent estimates. 

Table [5] lists the measured AB magnitudes (corrected 
to approximate total magnitudes through an aperture cor- 
rection) for the IRAC 3.6 /im, 4.5 /jm and 5.6 /im channels. 
None of the sources were detected at > 3 a in the 5.6 & 8 fim 
channels, although SBM03#3 has a very marginal (2 a) 
detection at 5.6 fim (AB = 24.4 ± 0.7). In the short- 
wavelength channels, SBM03#1 & #3 are well-detected (> 
10(7, AB m 24), and we have a more marginal 3 a detec- 
tion of GLARE#3001 (z' AB = 26.4) at AB » 26. Only 
SBM03#1 is detected at 4.5 ^m - SBM03#3 falls outside 
the region surveyed so far in this filter, and GLARE#3001 
is undetected. The fainter GLARE#3011 (z' AB = 27.2) was 
undetected at all IRAC wavelengths. For the HST photom- 
etry, we use the ACS i'-band and z'-band magnitudes from 
Stanway, Bunker & McMahon (2003) from GOODS-South 
and, in the case of SBM03#1, from the deeper UDF (Bunker 
et al. 2004). These magnitudes are already corrected to 
a total flux, also through an aperture correction. We use 
NICMOS magnitudes (Thompson et al. 2005) in F110W 
and F160W ('J' and 'H' band) from Stanway, McMahon 
& Bunker (2005) for SBM03#1, and for the others we used 
vl.O of the ESO VLT/ISAAC GOODS/EIS images 6 in the 
J and K s bands (Vandame et al. in prep). In the ground- 
based near-infrared ISAAC images, we used l"-diameter 
apertures. The seeing varied across the ISAAC field, as 
different tiles were taken over many nights, so we deter- 
mined the aperture correction from unresolved sources in 
each tile. For the J- and _ff s -band images the seeing is typi- 
cally good (FWHM = 0'.'4 - 0'.'5), and the aperture correc- 
tion is « 0.3 — 0.5 mag, determined from bright but unsatu- 
rated isolated stars measured in 6"-diameter apertures. 

The galaxy SBM03#3 is towards the edge of the re- 
leased ESO/VLT imaging, where fewer frames overlap and 
the noise is higher. Its K s magnitude seems to be anoma- 
lously faint (a 3 a detection at Kab = 25.5, twice as faint 
as the shorter- wavelength z' and J, which have AB = 
24.7 mag, and four times fainter than the longer-wavelength 
IRAC 3.6 ^m detection at AB = 24.0). We still include the 
K s filter in the fitting of the stellar populations fSection l2.3H . 



B Available from 

http : //www. eso . org/science/goods/releases/20040430/ 



but caution that the K s magnitude may not be reliable. For- 
tunately, this filter makes minimal difference to the best-fit 
populations as its statistical weighting is low. 

2.3 Spectral Energy Distributions 

The final step in the reduction process is the construction of 
spectral energy distributions (SEDs) for the chosen sources. 
Spitzer photometry and SED fitting has already produce 
great insight into the z = 2 — 3 Lyman-break star-forming 
galaxies (Barmby et al. 2004; Shapley et al. 2005) and here 
we apply similar techniques at z ~ 6. In order to compare 
our photometry with stellar spectral synthesis models, we 
utilise the latest Bruzual & Chariot (2003, hereafter B&C) 
isochrone synthesis code. We use the Padova evolutionary 
tracks (preferred by B&C 2003). The models utilise 221 age 
steps from 10 5 to 2 x 10 10 yr, approximately logarithmically 
spaced. Models with Salpeter (1955) initial mass functions 
(IMF) were selected, although in Section 13.51 we also con- 
sider the effect of adopting a Chabrier (2003) IMF. There 
are 6900 wavelength steps, with high resolution (FWHM 
3 A) and lA pixels over the wavelength range of 3300 A 
to 9500 A (unevenly spaced outside this range). From the 
full range of metallicities offered by the code, we considered 
both solar and l/5th solar models. From several star for- 
mation histories available, a single stellar population (SSP 
- an instantaneous burst), a constant star formation rate 
(SFR), and exponentially decaying (t) SFR models were 
used. This latest version of the Bruzual & Chariot (2003) 
models includes an observationally-motivated prescription 
for thermally-pulsing asymptotic giant branch (TP-AGB) 
stars (Liu, Graham & Chariot 2002), and also the TP-AGB 
multi-metallicity models of Vassiliadis & Wood (1993). 

For each of the four i'-drops with spectroscopic red- 
shifts, the filters were corrected to their rest-frame wave- 
lengths by the appropriate redshift factor. The measured 
flux was folded through the filter transmission profiles, and 
the best-fit age model was computed by minimising the re- 
duced x 2 ) using the measured errors on the magnitudes. 
The number of degrees of freedom is the number of inde- 
pendent data points (magnitudes in different wavebands). 
The Bruzual & Chariot spectra are normalized to an initial 
mass of 1 Mq for the instantaneous burst SSP model, and an 
SFR of 1 Mq yr -1 for the continuous star formation model. 
The fitting routine returned the normalisation for the model 
which was the best-fit to the broad band photometry (i.e., 
minimized the reduced X 2 ) - this normalisation was then 
used to calculate the corresponding best-fit stellar mass (see 
Section I3.3H . 

For SBM03#T, the SED-fitting process was conducted 
using photometry from the detections in the z' , F110W, 
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Figure 2. Images of SBM03#1 (z = 5.83), taken with HST/ ACS 
and z'-band); HST/NICMOS (1.1 (im F110W 'J-band', and 
1.6 /Ltm F160W '//-band'), VLT/ISAAC iY s -band (smoothed with 
a 3-pixel 0'.'45 boxcar) and the 4 Spitzer channels (3.6— 8.0 fim). 
Each panel is 8 arcsec across (a projected distance of 50kpc). 




Figure 3. Images of SBM03#3 {z = 5.78), taken with HST/ ACS 
(«'- and z'-band); VLT/ISAAC J- and K s -bands and Spitzer 
channels 3.6 & 5.6 /im. Each panel is 8 arcsec across (a projected 
distance of 50kpc). The J and K B band images have been 
smoothed using a 3-pixel box. 



F160W, K s , 3.6 /jm & 4.5 /xm filters, whilst for SBM03#3 
the z' , J, Ks, 3.6 /xm and 5.6 /xm data were used. 

Although some of our data points (particularly from 
the HST j ACS imaging) have S/N > 10, we set the mini- 
mum magnitude error to be A(mag) = 0.1 to account for 
flux calibration uncertainties. During the SED-fitting pro- 
cess, the i'-band flux was ignored, as this band is prone to 
contamination due to Lyman-a forest absorption shortwards 
of Lyman-a (A rcst = 1216 A) and also emission line contam- 
ination due to Lyman-a itself. 



Figure 4. Spitzer channel 1 (3.6 /im) image of SBM03#3 (right, 
lower-left of panel), with the best-fit GALFIT model for the sur- 
face brightness profile of the nearby neighbour (middle image, an 
exponential disk convolved with the PSF), and the Spitzer image 
with the neighbouring galaxy subtracted (left). The photometry 
of SBM03#3 is only affected at the ~ 5 per cent level by confusion 
from the neighbour. The frame is 9 arcsec across. 




Figure 5. Images of GLARE#3001 (z = 5.79), taken with 
HST/ ACS and z'-band); VLT/ISAAC J- and 7\ s -bands and 
Spitzer channels 3.6 & 4.5 fim. Each panel is 8 arcsec across (a pro- 
jected distance of 50kpc). The J and K s band images have been 
smoothed using a 5-pixel box. 



3 ANALYSIS 

For two of the spectroscopically-confirmed galaxies 
(SBM03#1 & #3) we have robust Spitzer detections at 
3.6 /im of AB « 24mag (> 10 <r), with the fainter GLARE 
sources marginally detected (GLARE#3001) or undetected 
(GLARE#3011). Progressing in increasing wavelength from 
the HST/ ACS z'-band and near-infrared (0.9 - 2.2 fim) to 
the Spitzer/IRAC 3.6 fim channel, the two well-detected 
i'-drops brighten by A(mag)AB = 0.8 (Figures HJ&Q - a 
factor of 2 in flux density, /„. SBM03#1 is also robustly 
detected with a similar magnitude at 4.5 fim. We now 
consider the implications of the spectral energy distribu- 
tions of these i'-drops, specifically which star formation 
histories can produce the observed spectral breaks between 
2.2 - 3.6 fim. 



3.1 Balmer/4000 A Breaks in z w 6 Galaxies 

The presence of a Balmer/4000 A break is suggestive of a 
system viewed a significant time after a major epoch of star 
formation - the Balmer break at 3648 A is strongest in an A- 
star population, and metal line blanketting (predominantly 
Fe n) due to an older late-type stellar population produces 
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ID 


i' 


z' J K s 3.6 fim 


4.5 fim 


5.6 fim 




SBM03#1 


26.99 ±0.03 


25.35 ±0.02 25.54t±0.04 25.14 ±0.22 24.24 ± 0.09 


24.42 ±0.15 


> 23.8 (3 a) 




SBM03#3 


26.27 ±0.13 


24.67 ±0.03 24.72 ±0.15 25.55 ± 0.40* 23.94 ± 0.07 




24.39 ±0.69 




GLARE#3001 28.03 ±0.14 


26.37 ±0.06 26.11 ±0.32 25.22 ± 0.32 25.88 ± 0.36 


> 26.1 (3 a) 


> 23.8 (3a) 




GLARE#3011 > 28.8 (3 a) 


27.15 ±0.12 >26.2(3cr) >25.6(3cr) >26.5(3cr) 


> 26.1 (3cr) 


> 23.8 (3 a) 






° The K s 


magnitude of SBM03#3 is anomalously faint. 








t From NICMOS F110W imaging (Stanway, McMahon & Bunker 2004), rather than ESQ VLT/ISAAC imaging. The NICMOS F160W 








magnitude of SBM03#1 is H A B = 25.51 ± 0.05 








Table 2. Magnitudes (AB system) of four i'-band drop galaxies. SBM03#3 is not in the field of view for the 4.5 fim 


and 8.0 /im filters. 




All four galaxies 


are undetected at f 


i.O/im, with AB > 23.5 (3 a). 








Model 


Reduced xLn 


Age / Myr E(B-V) M tot / 1O 1() M M stollar / 10 10 M Q 


Current SFR / Mq yr" 1 b / 10 


-io yr -i 


burst (SSP) 


2.83 


90.5 0.00 1.35 1.13 


0.00+ 





.00 


r = 10 Myr 


2.81 


102 0.00 1.34 1.14 


0.05+ 





.04 


t — 30 Myr 


2.34 


143 0.02 1.63 1.36 


4.59 


3 


.38 


t — TO Myr 


1.61 


255 0.00 2.05 1.68 


7.88 


4 


69 


t = 100 Myr 


1.38 


321 0.00 2.33 1.97 


9.80 


4 


.20 


T = 300 Myr 


0.94 


641 0.00 3.41 2.70 


15.2 


4 


.97 


T = 500 Myr 


0.90 


905° 0.00 4.35 3.38 


17.0 


5 


.03 


t = 1000 Myr 


0.97 


1280° 0.00 5.46 4.21 


21.1 


5 


.01 


continuous 


1.15 


5000° 0.00 7.93 5.84 


15.9 


2 


.72 




t Ruled out because SFR below lower limit set by Lyman-a emission. 








Ruled out as age is close to or exceeds 1 Gyr, age of Universe at z : 


5i 6 






Table 3. Favoured model parameters of various SED fits for SBM03#1, mctallicity = Zq. We varied age between and 20 Gyr, and 




extinction between 0.00 and 1.00. We list both the mass currently in stars (M sto n ar ), and the total baryonic mass (M tota j which includes 




the mass returned to the IGM by evolved stars) for each star formation history. 








Model 


Reduced X 2 min 


Age / Myr E(B-V) M tot / 1O 1(, M M stollar / 1O 1O M 


Current SFR / M yr" 1 b / 10 


-io yr -i 


burst (SSP) 


1.54 


114 0.00 1.32 1.09 


0.00+ 





.00 


r = 10 Myr 


1.50 


128 0.00 1.35 1.13 


0.02+ 





.02 


r = 30 Myr 


1.30 


161 0.01 1.42 1.17 


2.23+ 


1 


.91 


t = 70 Myr 


0.99 


255 0.00 1.59 1.29 


6.12 


4 


.74 


t = 100 Myr 


0.87 


321 0.00 1.83 1.47 


7.70 


5 


.24 


r = 300 Myr 


0.77 


641 0.00 2.77 2.18 


12.4 


5 


.69 


t = 500 Myr 


0.78 


905° 0.00 3.57 2.76 


14.0 


5 


.07 


t = 1000 Myr 


0.81 


1280° 0.00 4.54 3.48 


17.5 


5 


.03 


continuous 


0.93 


2600° 0.02 4.34 3.35 


16.7 


4 


.99 




t Ruled out because SFR below lower limit set by Lyman-a emission. 








Ruled out as age is close to or exceeds 1 Gyr, age of Universe at z : 


% 6 






Table 4. Favoured model parameters of various SED fits for SBM03#1, mctallicity = O.2Z . We varied age between and 20 Gyr, and 




extinction between 0.00 and 1.00. We list both the mass currently in stars (M ste n ar ), and the total baryonic mass (M tota j which includes 




the mass returned to the IGM by evolved stars) for each star formation history. 








Model 


Reduced X 2 min 


Age / Myr E(B-V) M tot / 1O 1(, M M stollar / 1O 1O M 


Current SFR / M yr" 1 b / 10 


-lOyr-i 



burst 3.66 47.5 0.00 1.12 0.97 0.00 f 0.00 

r = 10Myr 3.65 64.1 0.00 1.17 1.02 1.93+ 1.89 

r = 30Myr 3.21 114 0.00 1.53 1.29 11.7 9.07 

r = 70Myr 2.75 203 0.00 2.22 1.84 18.6 10.1 

T= 100 Myr 2.55 255 0.00 2.54 2.08 21.5 10.3 

t = 300 Myr 2.06 453 0.00 3.54 2.85 33.4 11.7 

T = 500 Myr 1.96 641 0.00 4.78 3.77 36.8 9.76 

t = 1000 Myr 1.88 905° 0.00 6.81 5.37 46.3 8.62 

continuous 1.82 1800° 0.00 5.20 4.02 28.9 7.19 

t Ruled out because SFR below lower limit set by Lyman-a emission. 
Ruled out as age is close to or exceeds 1 Gyr, age of Universe at z Rj 6 



Table 5. Favoured model parameters of various SED fits for SBM03#3, mctallicity = Zq. We varied age between and 20 Gyr, and 
extinction between 0.00 and 1.00. We list both the mass currently in stars (M ste n ar ), and the total baryonic mass (M tota j which includes 
the mass returned to the IGM by evolved stars) for each star formation history. 
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Model 


Reduced xL^ 


Age / Myr 


E(B-V) 


A/tot / 1O 1O M 


A^stollar / 10 10 M O 


Current SFR / A^yr" 1 


b / lO-^yr" 1 


burst 


2.15 


71.9 


0.00 


1.38 


1.17 


o.oot 


0.00 


t — 10 Myr 


2.12 


80.6 


0.00 


1.30 


1.12 


0.4lt 


0.37 


t — 30 Myr 


1.92 


128 


0.00 


1.51 


1.26 


7.21 


5.72 


t = 70 Myr 


1.73 


203 


0.00 


1.81 


1.50 


15.2 


10.1 


t = 100 Myr 


1.67 


255 


0.00 


2.07 


1.69 


17.5 


10.4 


T = 300 Myr 


1.50 


453 


0.00 


2.96 


2.37 


27.9 


11.8 


t = 500 Myr 


1.46 


571 


0.00 


3.65 


2.90 


34.2 


11.8 


t = 1000 Myr 


1.41 


806° 


0.00 


5.37 


4.22 


43.3 


10.3 


continuous 


1.36 


1400° 


0.00 


4.34 


3.44 


25.1 


7.30 



t Ruled out because SFR below lower limit set by Lyman-Q emission. 
Ruled out as age is close to or exceeds 1 Gyr, age of Universe at z ss 6 

Table 6. Favoured model parameters of various SED fits for SBM03#3, metallicity = 0.2Zq. We varied age between and 20 Gyr, and 
extinction between 0.00 and 1.00. We list both the mass currently in stars (M ste n ar ), and the total baryonic mass (M tota i which includes 
the mass returned to the IGM by evolved stars) for each star formation history. 
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Figure 6. Best-fit Bruzual & Chariot model for SBM03#1: an 
exponentially decaying star formation rate with r = 300 Myr, 
viewed 640 Myr after the onset of star formation. The stellar mass 
is 3.4 X 10 10 Mq . Flux density is in /„ units. 



Figure 7. Best-fit Bruzual & Chariot model for SBM03#3: an 
exponentially decaying star formation rate with r = 500 Myr, 
viewed 640 Myr after the onset of star formation, and with a stel- 
lar mass of 4.8 X 10 10 Mq. Flux density is in /„ units. 



the 4000 A break. For the two i'-drops with the best Spitzer 
detections (SBM03#1 & #3) we find evidence of a signifi- 
cant Balmer/4000 A break: the brightening by a factor of 2 
in /„ from the near-infrared (« 0.9 — 2.2 fim) to 3.6 /im im- 
plies a break amplitude of ~ 1.7 (for f\ flux densities). The 
z' , J&zK 3 colours are relatively flat in /„ (Figures 
ignoring the discrepantly faint K 3 magnitude of SBM03#3- 
which favours the spectral break interpretation rather than 
dust reddening (see Section 13.51 . The break amplitude of 
~ 1.7 is comparable to that seen by Le Borgne et al. (2005) 
in massive post-starburst galaxies at much lower redshifts 
(z ~ 1) from the "Gemini Deep Deep Survey" project. In- 
deed, the break amplitude in our two significant z ~ 6 cases 
is only slightly less than the D4000 (Bruzual 1983) observed 
at z ~ in the Sloan Digital Sky Survey (e.g., Kauffmann 
et al. 2003), with D4000 ft! 1.7 - 2.0. 

For all the galaxies discussed here, there has been re- 



cent or ongoing formation of at least some massive stars 
prior to the epoch of observation, as our Keck/Gemini spec- 
tra show Lyman-a emission produced through photoioniza- 
tion of hydrogen by the short-lived OB stars (Bunker et 
al. 2003; Stanway et al. 2004a). In order to produce the 
Balmer/4000 A break amplitude at A rcst ~ 4000 A, most of 
the stellar mass probably formed well before the current 
starburst, most likely > 100 Myr previously f Section 13.21 . 
Hence the galaxies SBM03#1 & #3 are already, to some ex- 
tent, established systems. This is a significant finding since, 
at z ~ 6, the Universe is less than 1 Gyr old. 

There is no significant evidence for a Balmer/4000 A 
break in the two (fainter) spectroscopically-confirmed i! - 
drops: the 3 a detection of GLARE#3001 at 3.6 (im sug- 
gests a break of ± 0.5 mag, and the non-detection of 
GLARE#3011 places a 3 a upper limit of AAB = 1 mag on 
the break amplitude. Hence the constraints for the stellar 
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0.6 F 




Figure 8. Filter profiles of the four Spitzer infrared wavebands. 




1 2 3 4 5 

Figure 9. The best-fitting two-component stellar population 
model (Salpeter IMF) for SBM03#1: a dominant 450 Myr popu- 
lation of mass 3.6 X 10 10 Mq, with some ongoing star formation 
activity (a burst for the last 10 Myr involving 0.7% of the stellar 
mass). Using a Chabrier IMF produces an identical best-fit age, 
with a similar mass fraction in the burst (0.6%), but a 30% lower 
total stellar mass (2.5 X 10 10 Mq). 

ages of GLARE#3001 & 3011 are weak, but they are consis- 
tent with younger stellar populations than SBM03#1 & #3, 
which are brighter and presumably more massive. In the re- 
mainder of this discussion, we focus on these two brighter 
sources which have significant Spitzer detections for which 
we can estimate the stellar ages and masses. 

3.2 Ages of the i'-drop galaxies 

We explored the best-fit stellar ages and star formation 
histories of the galaxies SBM03#1 & #3 by comparing the 
SEDs from our broad-band photometry with the popula- 



tion synthesis models of Bruzual & Chariot (2003). We 
cover a range of star formation histories (presented in Ta- 
blcs|3|EIinilSI&0 designed to bracket most plausible evolu- 
tion scenarios. Only a subset of these star formation histories 
provided acceptable fits to our photometry. We now discuss 
several classes of model and their validity. 

We begin by considering an idealized "simple stellar 
population" (SSP) model, where a galaxy is viewed some 
time after its entire stellar mass formed in an instantaneous 
burst. These models provide relatively poor fits (reduced 
X 2 ~ 3 for solar metallicity) to our data. This is under- 
standable given that we know that there is at least some 
ongoing star formation in both galaxies. The SSP models 
provide an absolute lower age limit for the bulk of the stel- 
lar mass. These SSP models yield an age of » 100 Myr for 
the best fit population (and > 70 Myr at 99% confidence) in 
the absence of dust reddening (see Section f3.5|l . implying a 
formation redshift of Zf ^> 6.4. 

At the other extreme, a constant star formation history 
was considered. The best fit ages in this simplistic scenario 
are 1.5-5 Gyr for metallicity 0.2-1 Z @ (with SBM03#3 best- 
fit with a younger population than SBM03#1, and lower 
metallicities decreasing the age); these ages are obviously 
unphysical, as they exceed the age of the universe at this 
redshift (« 1 Gyr). This strongly implies that we are seeing 
both galaxies at an epoch when the star formation rate is 
declining, or that the current burst was preceded by a more 
significant episode of star formation. 

Hence we next considered a range of more realistic mod- 
els including an exponentially decaying star formation rate 
(SFR), such that the current SFR t = SFR e" t/T , where 
SFRo is the star formation rate at the onset of the burst and 
t is the decay time in Myr. Such models are intermediate be- 
tween SSP (r — * 0) and continuous star formation (r — > oo) 
models. We considered several values of r ranging between 
10 Myr to 1 Gyr. As before, those models with ages older 
than the Universe at z » 6 were rejected. Also, models with 
t — 10 Myr were disregarded as they yielded current SFRs 
below the lower limit set by the Lyman-a emission (Sec- 
tion 13. 6t . The favoured models have an age of » 600 Myr, 
with decay timescales of r = 300 Myr for SBM03#1, and 
r = 500 Myr for SBM03#3. 

Finally, we consider a model where the galaxies are com- 
posed of two distinct stellar components: an ongoing star- 
burst at the time of observation and an older population 
which formed in an instantaneous burst some time previ- 
ously. We considered starbursts with a constant SFR which 
started 3, 10, 30, & 100 Myr prior to the epoch of observa- 
tion. We varied the ratio of total stellar masses in these two 
populations, and found the best-fit ages for the old compo- 
nent to be 400-500 Myr. From the best-fit models (Table0, 
the fraction of total stellar mass being formed in the star- 
burst is 0.5 — 5% for starbursts of 3 — 100 Myr duration. 

In summary, our SED fitting shows that the broad-band 
colours can be fit with a variety of stellar ages/star formation 
histories. The lower limit on the age is > 100 Myr (from 
an SSP model) with the oldest allowed models comparable 
to the age of the Universe at z ~ 6. Our best-fits to the 
broad-band photometry come from exponentially-decaying 
star formation histories, or a two-component model where 
only 0.5 — 5% of the stellar mass is forming in an ongoing 
starburst. These best-fit models have mean stellar ages of 
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260 — 640 Myr, which would require formation redshifts of 
z f » 7.5 - 13.5. 

It has been speculated that certain treatments of the 
thermally-pulsing asymptotic giant branch (TP-AGB) may 
lead to different age estimates from broad-band colours. 
Maraston (2005) suggests the z ~ 2 — 3 population of 
'iEROs' identified by Yan et al. (2005) in HST+Spitzer 
images can be fit by younger ages than indicated by the 
Bruzual & Chariot (2003) models. Potentially the TP-AGB 
has a significant impact on the restframe V — K colours mea- 
sured in the ACS/IRAC images at z ~ 2 — 3. However, our 
i'-drops are at higher redshifts, and the TP-AGB has little 
effect at A < 0.7 /im (see fig. 14 of Maraston 2005), the wave- 
length range which we are probing with our HST+Spitzer 
detections. 



3.3 Stellar Masses at z « 6 

The SED fitting procedure described above leaves the nor- 
malisation of the Bruzual & Chariot model as a free pa- 
rameter (along with the age for a particular input model). 
Our code outputs this normalization and, by using the lu- 
minosity distance to the i'-drop galaxies, we can calculate 
the corresponding best fit stellar mass for each star forma- 
tion history. Tables |SJ HI El & 111 give the best-fit masses re- 
turned by the SED fitting. We list both the mass currently 
in stars (M ate n ar , and the total baryonic mass (M tota i which 
includes the mass returned to the IGM by evolved stars) 
for each star formation history. To assess the errors on the 
stellar mass estimates, we took the best-fit mass, M s t e iiar, 
and recalculated the reduced \ 2 f° r the same age, metal- 
licity and star formation history, but using total masses in 
the range 0.1 — 3M ste n ar (shown in Figure ITO1 . In all our 
models, the stellar masses of SBM03#1 & 3 were > 10 10 M Q 
at 95% confidence (2 a). The lowest masses were returned 
by the SSP model, and the exponentially-decaying models 
with short decay times (r ~ 10 Myr): in fact, these models 
are unable to provide sufficient ongoing star formation to 
explain the lower limit set on the star formation rate by the 
Lyman-Q emission (see Section I3.6H . Our preferred models 
(r = 70 — 500 Myr and the two-component stellar popula- 
tion) have stellar masses of « 2 - 4 X 10 10 M Q for a Salpeter 
IMF (we consider the impact of adopting a Chabrier IMF 
in Section 1331 . 

We measure stellar masses from the best-fitting SEDs 
of M > 2 x 10 10 Mq at z « 6. This is surprisingly large, 
supporting our contention that at least these two objects 
are well-established galaxies. The stellar mass is equivalent 
to 20% of that for a L* galaxy today, using L* — —21.21 
from the SDSS analysis of Blanton et al. (2003) and tak- 
ing M/Ly ~ 5Mq/Lq (appropriate for a ~ lOGyr old 
population from B&C models using a Salpeter IMF) to ob- 
tain M* — 1.2 x 1O 11 M0, comparable to the estimate of 
M* = 1.4 x 10 11 M Q from Cole et al. (2001) for our adopted 
Salpeter IMF. 



3.4 Stellar Mass Density and the Kormendy 
Relation 

A key question given the increased information content on 
our two z ~ 6 sources is the likely descendant population 



and connection with systems seen at lower redshift. In this 
respect, we now consider the high surface brightness and 
small angular extent of the i'-drop galaxies. 

To facilitate this discussion, we examine the properties 
of the z k 6 galaxies with reference to the z ~ Kor- 
mendy relation between galaxy surface brightness and half- 
light radius (Kormendy 1977). Both SBM03#1&#3 have 
half-light radii (effective radii) of Vhi ~ 0'.'08 (Bunker et 
al. 2004; Bunker et al. 2003), and IRAC total magnitudes 
of AB w 24 at 3.6 /im (corresponding to the rest-frame V- 
band). The average surface brightness within the effective 
radius, (I)hl, is related to the average surface brightness at 
ruu hi, by (hi) = 3.6 x (I) M for a de Vaucouleurs profile 
(Ciotti & Bertin 1999). Hence our surface brightness at ru 
is /j, e — 21.6 mag arcsec -2 . 

Cosmological surface brightness dimming of (1 + z) 4 
would have dimmed real surface brightnesses by +8.3 mag, 
and passive evolution from an age of w 200 Myr at z — 5.8 
to 12.8 Gyr at z — would produce a luminosity dim- 
ming of a factor of x30 in V (+3.7 mag, using the B&C 
models). Hence, we might expect SBM03#1&#3 to have 
surface brightnesses today of/ic = 21.6 — 8.3 + 3.7 = 
17.0 mag arcsec -2 for their half-light radii of ru = 0.4 kpc. 
In fact, the passively-evolved surface brightnesses sit com- 
fortably on the present-day Kormendy relation (e.g., that 
of Ziegler et al. 1999 for cluster ellipticals), extrapolated to 
smaller scale lengths. Hence it seems that the inferred prop- 
erties of these z ~ 6 galaxies are compatible with galaxy 
scaling relations at z = 0, subject to dimming through stel- 
lar evolution. Although our analysis has focused on only 
two of the most luminous systems, chosen by virtue of their 
luminosity, the possible implications are profound. What 
do these i'-drop galaxies at z ~ 6 evolve into? Given that 
these objects are barely resolved in the HST /ACS data with 
rhi < 0.5 kpc, they are unlikely to passively evolve into the 
ellipticals we see in the present-day Universe - these are 
typically larger. Merging would be required to explain the 
size evolution, although we note that the stellar masses and 
spatial scales of our z w 6 galaxies are similar to those of 
some spiral bulges today: a stellar age of w 12 Gyr (Wyse, 
Gilmore & Franz 1997) today would imply a formation epoch 
of z > 5. 

3.5 Effects of Metallicity, Dust and IMF on SED 
Fitting 

When conducting the SED fitting, two metallicities were 
considered: solar (Zq) and a sub-solar model (0.2.^0 ). The 
ages and masses of the best-fit models were similar for both 
metallicities, although the sub-solar model returned slightly 
better fits to the data, with smaller reduced Xmin values. 

We also considered whether the red optical-infrared 
colours (spanning the rest-frame UV to optical) could be 
attributable to dust reddening instead of, or as well as, an 
underlying old stellar population. To further examine this, 
we adopt the reddening model of Calzetti (1997) 7 , appropri- 

7 The Calzetti reddening is an empirical law is given in terms of 
the colour excess E(B — V) = 0.44_B gas (B — V) with the wave- 
length dependence of the reddening expressed as: 

k(X) = 2.656(-2.156 + 1.509/A - 0.198/A 2 + 0.011/A 3 ) + 4.88 
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burst duration / Myr reduced \ 2 f ra c of mass current SFR / Mq yr 



age / Myr stellar mass / 10 10 Mq 



3 
10 
30 
100 



0.71 
0.68 
0.72 
0.79 



0.5% 
0.7% 
1.7% 
5.2% 



60.0 
25.3 
19.0 
68.3 



400 
450 
450 
450 



3.6 
3.6 
3.3 
3.1 



Table 7. Two-population composite models for SBM03#1, assuming an ongoing burst of constant star formation rate, commencing 
3-100 Myr ago, within an older galaxy. The best- fit galaxy age and mass fraction of the starburst for each model are tabulated (for solar 
mctallicity, Zq). 




age / yr age / yr 



Figure 10. The allowed range of masses for several exponentially- 
decaying SFR, models for SBM03#1, with decay times ranging 
from t — 10 — 1000 Myr, as well as an instantaneous burst model 
and a constant SFR model. Those with stellar ages > 1 Gyr (right 
shaded region) are ruled out (the Universe is younger than this at 
z ~ 6). The models with ages < 10 s yr are also excluded as they 
fail to produce the observed current SFR inferred from Lyman-a 
emission (left shaded region). Contours are 68% confidence (solid 
line), 95% confidence (dashed line) and 99% confidence (dotted 
line) for reduced \ 2 of ss 1, 2, 3. 

ate to starburst galaxies. We constructed a grid of models for 
both the SSP (instantaneous burst) and constant SFR sce- 
narios, using the Bruzual & Chariot templates over the full 
range of 221 age steps. For each age, we varied reddening in 
the range E(B - V) = 0.00 - 1.00 mag, in steps of 0.01 mag, 
computing the reduced \ 2 a t each step. We find little evi- 
dence for substantial dust reddening in the detected starlight 
out to « 5 fim: the formal best-fits for both SBM03#1 & #3 
are no reddening. This is consistent with the flat spectral 
slopes in /„ between the ACS z'-band and the NICMOS 

for 0.12 fj,m < A < 0.63 /im, and 

fc(A) = [(1.86 - 0.48/A)/A - 0.1]/A + 1.73 

for 0.63 /an < A < 1.0 fim, and the flux attenuation is given by: 



Figure 11. SBM03#1: A plot showing the reddening E(B-V) 
values for the two limiting scenarios of an instantaneous burst 
model and a constant SFR, model, for the case metallicity = Zq. 
Contours are 68% confidence (solid line), 95% confidence (dashed 
line) and 99% confidence (dotted line) for reduced x 2 °f ~ 1> 2, 3. 



J- and //-bands reported by Stanway, McMahon & Bunker 
(2005) for many of the i'-drops in the Hubble Ultra Deep 
Field, including SBM03#1. The brightening in f v flux be- 
tween the near-infrared (0.8 — 2.2 (im) and the IRAC bands 
at 3 — 5 fim is best explained by a spectral break rather than 
the smoother continuum gradient produced by dust redden- 
ing, and this is reflected in the best-fit stellar populations. 

Finally, we tested the effects on the derived stellar 
masses of the assumed initial mass function. For the galaxy 
SBM03#1 (where the multi- waveband data was best, as 
it included NICMOS imaging in the near- infrared), we re- 
ran the 'two-population' model fits (Section 13.21 with a 
Chabrier (2003) IMF, instead of a Salpeter (1955) power 
law IMF. We used the same ongoing 10 Myr constant star 
formation rate burst, and fit for the age of an underlying 
older stellar population and the relative stellar masses of 
the burst and the old stars. Both IMFs produced compara- 
bly good fits (near- identical % 2 values), the same 450 Myr 
ages, and near-identical burst fractions (0.6% and 0.7% by 
mass for the Chabrier and Salpeter IMFs). The main differ- 
ence came in the best-fitting total stellar mass: the Chabrier 
model produced a mass « 30% less than the Salpeter IMF 
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Figure 12. As Figure [HI but for SBM03#3. Contours are for 
reduced \ 2 of ~ 1>2,3, 4, corresponding to 68% confidence (solid 
line), 95% confidence (dashed line), 99% confidence (dotted line), 
and 4 a confidence (dot-dash line) . 

(2.5 x 10 10 M Q compared with 3.6 x 10 10 M Q for SBM03#1). 
This effect is mainly a mass re-scaling independent of star 
formation history, as the discrepancy primarily arises from 
the different mass fraction in long-lived low-mass stars. 



3.6 Star Formation Rates 

We can usefully compare the ongoing star formation rate, 
the stellar mass and implied age, to deduce whether our two 
selected z m 6 galaxies are being seen during a fairly quiet 
or active period in their history. 

First, let us consider the ongoing SFR. Star forma- 
tion will dominate the rest-frame UV light (probed by the 
HST/ ACS images) in the absence of dust obscuration or 
a significant AGN contribution. In our first analyses of £'- 
drop galaxies in the GOODS ACS images (Stanway, Bunker 
& McMahon 2003; Stanway et al. 2004a) and the Ultra Deep 
Field (Bunker et al. 2004) , we used the conversion from rest- 
frame UV flux density to SFR, given by Madau, Pozzetti & 
Dickinson (1998), appropriate for continuous star formation. 
Based on the a'-band magnitudes we inferred unobscured 
SFRs of 19.5 & 33.8 Meyr -1 for SBM03#1&#3, after ac- 
counting for Lyman forest blanketting of Da ~ 0.95 short- 
wards of Lyman- a. 

The fits of the Bruzual & Chariot stellar synthesis mod- 
els to the broad-band photometry provide estimates of the 
current SFR for a range of histories (TablesEHUEIEI&Q . 
However, a firm lower limit arises from our Keck/DEIMOS 
spectroscopy (Bunker et al. 2003, Stanway et al. 2004a), 
which revealed Lyman-a emission. If this line emission is 
powered by ionizing photons from OB stars, then there must 
be star formation activity in these galaxies within the past 
10 Myr (the lifetime of these stars, which dominate the Ly- 
man continuum flux). 



We have argued previously that the absence of X-ray 
emission and high-ionization lines such as Nv 1240 A , cou- 
pled with the narrow velocity width of the Lyman-a emis- 
sion, renders an AGN interpretation unlikely. Our spec- 
troscopy of SBM03#1&#3 indicates line fluxes of 2 x 
10~ 17 erg s -1 cm -2 for both sources, with rest-frame equiv- 
alent widths of W Te st ~ 30, 20 A. Assuming case B recom- 
bination and the same Salpeter IMF as in the Bruzual & 
Chariot models, the star formation rates from Lyman-a are 
ft* 6 Mq yr _1 . This may be treated as a firm lower limit on 
the current SFR, as the resonantly-scattered Lyman-a line 
is invariably quenched by dust to well below its case B line 
strength. 

Our spectroscopic lower limit of an SFR > 6M©yr _1 
rules out the simple SSP model and those declining star 
formation rate models with the shortest decay times (r = 
10 — 30 Myr, Figure I1UI . Our favoured models are an expo- 
nentially decaying SFR with decay time r ~ 300 — 500 Myr, 
and two-component models with 1 — 2% of the stellar mass 
created in an ongoing starburst which began 10 — 100 Myr 
prior to the epoch of observation. These all indicate similar 
SFRs to our original estimates of « 20 — 30 Mq yr _1 from 
the z'-band flux. 

Returning to the ratio of the current SFR to the stel- 
lar mass already formed, it is helpful to consider the '6- 
parameter' (e.g., Brinchmann et al. 2004) - a measure of 
the fraction of the total stellar mass currently being born 
as stars, defined as b — SFR/Msteiiar- We measure b w 
5&10 x lO-^yr -1 for SBM03#1&#3, comparable with 
that inferred by Egami et al. (2005) from the Spitzer image of 
a lensed galaxy with a photometric redshift of z ~ 7 (Kneib 
et al. 2004), and an indication of vigorous current star for- 
mation which supports the conclusion of this paper. How- 
ever, the past-average SFR must actually be greater than 
the current SFR in order to build our « 2 x 10 10 M Q galax- 
ies, given the short time available prior to z ~ 6 (lGyr). 
This is why decaying SFR models appear to give the best 
fits to the SEDs. Juneau et al. (2005) define a character- 
istic growth timescale (in a study of z ~ 1 galaxies), with 
tsFR = M s toiiar/SFR = 1/6. The ratio of this to the Hubble 
time at that epoch, i.sr(z), can be interpreted as galaxies in 
a declining or quiescent star formation mode if tsFR > tn- 
This is the case for both SBM03#1&3 for the best-fitting 
models. Indeed, to form the 2 — 4 x 10 10 Mq of stars re- 
quires an average star formation rate of 50 — lOOMQyr -1 
over the redshift range Zf = 7.5 — 13.5 favoured for the for- 
mation of the old stellar component (Section 13.21 . Hence, 
(SFR 3= 7.5_i3.5)/(SFR z=6 ) ~ 2 — 5; the past-average star 
formation rates of SBM03#1 & 3 are factors of 2 — 5 greater 
than their current rates at the epoch of observation [z ~ 
5.8). 

Although our chosen galaxies are the brightest con- 
firmed i'-drops in the GOODS-South field and possibly un- 
representative, the present work does imply that in these 
galaxies at least there was a yet earlier vigorous phase of 
activity, possibly at z > 10, which may have played a key 
role in reionizing the Universe. This may be consistent with 
the measurement of temperature-polarization correlation of 
the cosmic microwave background from the Wilkinson MAP 
satellite by Kogut et al. (2003). 
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4 CONCLUSIONS 

Our group previously identified and spectroscopically- 
confirmed z ~ 6 galaxies through HST/ ACS i'-drop imag- 
ing and Keck/DEIMOS & Gemini/GMOS spectroscopy. 
This paper presents the first infrared detections of this 
population using Spitzer/IRAG. We have significant (« 
10 a, AB ~ 24 mag) detections at 3.6 /im of the Stan- 
way, Bunker & McMahon (2003) galaxies #1 and #3 (at 
Zspcc = 5.83,5.78), and a more marginal detection of the 
•Zspcc = 5.79 galaxy GLARE#3001 (Stanway et al. 2004b). 
We also detect SBM03#1 at 4.5 ^m (SBM03#3 is outside 
the field of view of this 4.5 fim filter). 

We infer from Lyman-a emission in our discovery spec- 
tra that there is ongoing star formation of > 6MQyr _1 
(as would be expected in these rest-UV-selected objects). 
However, the preceding star formation history has not been 
explored until now. In the two best-detected galaxies, we 
have evidence of a significant Balmer/4000 A break, indica- 
tive of a prominent older stellar population which probably 
dominates the stellar mass. Exploring a range of popula- 
tion synthesis models indicates that the average stellar age 
is > 100 Myr; our best-fit models suggest preferred ages of 
250 — 650 Myr for an exponentially-declining star formation 
rate (of decay time r ~ 70 — 500 Myr) or a two-component 
model (with an ongoing starburst responsible for 0.5 — 5% 
of the total stellar mass). This implies formation epochs of 
z f w 7.5 - 13.5 for the galaxies SBM03#1 &#3. 

In all our models, the best-fit stellar masses are > 
10 10 M Q , with 95% confidence masses of 1.3-3.8 x 10 10 M Q . 
This indicates that at least some galaxies with stellar masses 
> 20% the mass of L* galaxies today were already assembled 
within the first Gyr of the Universe. For these objects, the 
past average star formation rate is comparable to, or greater 
than the current SFR, implying that there may have been 
even more vigorous episodes of star formation at higher red- 
shifts. These may have played a key role in reionizing the 
Universe, consistent with the earlier studies of Bunker et al. 
(2004) and Egami et al. (2005). 
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